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Abstract 

We derive semiclassical approximations for wavefunctions. Green's functions and expecta- 
tion values for classically chaotic quantum systems. Our method consists of applying singular 
and regular perturbations to quantum Hamiltonians. The wavefunctions. Green's functions 
and expectation values of the unperturbed Hamiltonian are expressed in terms of the spectral 
determinant of the perturbed Hamiltonian. Semiclassical resummation methods for spectral 
determinants are applied and yield approximations in terms of a finite number of classical tra- 
jectories. The final formulas have a simple form. In contrast to Poincare surface of section 
methods, the resummation is done in terms of the periods of the trajectories. 
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1 Introduction 



Semiclassical methods have a central role in the study of classically chaotic quantum systems, 
because they allow to explain properties of quantum systems on the basis of the underlying 
classical dynamics. A prominent example is the Gutzwiller trace formula that approximates the 
density of states by a sum over classical periodic orbits. Among its many applications are the 
explanation of periodic orbit signatures in quantum spectra and the semiclassical evaluation 
of spectral statistics [HIS]. For a semiclassical calculation of individual quantum energy levels, 
however, it is more efficient to use formulations in terms of spectral determinants instead of 
the density of states. Resummation methods for semiclassical spectral determinants allow to 
express contributions of long orbits in terms of short ones and yield approximations in terms 
of a finite number of orbits [31 HI El E]. As a consequence, a significantly smaller number of 
periodic orbits is needed. Numerical tests show that these resummed semiclassical sums can 
be remarkably accurate approximations to quantum spectral determinants [7]. 

Similar methods for wavefunctions and related quantities have been developed in several 
articles. Bogomolny derived a semiclassical theory for averages of wavefunctions over small 
windows of energy and position, starting from the semiclassical Green's function [^. A corre- 
sponding approach in phase space was developed by Berry One main motivation behind 
these theories was to explain the scar phenomenon [10] . They apply, however, to averages over 
several wavefunctions within the considered energy window and do not involve resummations. 
The first semiclassical formula for single wavefunctions was obtained by Agam and Fishman 
[TT] . Their formula is based on a semiclassical resummation and was used to accurately predict 
the occurrence of scars in wavefunctions In subsequent articles Fredholm methods for 

integral equations on quantum Poincare surfaces of section were applied to obtain simpler ex- 
pressions and put the work of Agam and Fishman on a more rigorous basis [T3l [HI [151 [HI [13 • 
In these approaches the resummation is done in terms of the number of intersections of the 
trajectories with the Poincare surface instead of their periods. In [TB] special properties of tiling 
billiards on the pseudo-sphere were used to obtain a resummation in terms of the periods for 
these systems. 

In this article we use a different approach. We apply perturbations to a Hamilton in or- 
der to express wavefunctions. Green's functions and expectation values in terms of spectral 
determinants. This allows the application of semiclassical resummation methods for spectral 
determinants. The resulting formulas are simpler than previous approximations, although they 
are equivalent to them in the semiclassical limit. 

The article is organized as follows. In section 2 we introduce spectral determinants and 
their semiclassical approximations. In section 3 we discuss singular rank one and rank two 
perturbations of Hamiltonians which are applied in section 4 to obtain resummed semiclassical 
approximations for wavefunctions and Green's functions. Corresponding results for expectation 
values are obtained in section 5, and section 6 contains the conclusions. 

2 Spectral determinants 

Spectral determinants are functions whose zeros are given by the eigenvalues of a quantum 
Hamiltonian. They are entire functions of the energy E, real valued for real E, and can be 
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expressed as a product over the quantum energy levels 

oo 

A{E) = II A{E, Em) {E - Em) . (1) 



m=l 



The coefficients A{E, Em) are non-vanishing functions of E and make the product convergent. 
The specific form of these coefficients depends on the chosen regularization. Two common 
forms of regularizing spectral determinants, the Weierstrass and the zeta regularization, are 
discussed in appendix El 

As is also discussed in this appendix, the semiclassical approximation to the spectral deter- 
minant can be obtained from the (regularized) trace of the Green's function and is expressed 
in terms of the trajectories of the corresponding classical system [19] 

A{E) ^ B{E)e-'^^^^\{E) . (2) 

Here B{E) is a function which is real for real E, see appendix |Al It will drop out of the final 
formulas of this article. N{E) is the mean spectral counting function which originates from 
zero-length trajectories. Its leading order asymptotic term is given by the Thomas-Fermi or 
Weyl estimate 

N{E) {2nh)-'^ J d'^qd^pOiE - H{q,p)) as ^ ^ , (3) 

where 0{E) denotes the Heaviside theta function and D is the dimension. 

The function C(-^) receives contributions from the classical periodic orbits. In the case of 
chaotic systems it can be expressed as an Euler product over the periodic orbits. Alternatively, 
the product can be expanded and expressed as a sum over composite (or pseudo) orbits. We 
discuss the second form in the following. To each composite orbit, labeled by n, corresponds 
a finite set of primitive periodic orbits P„, and the action of the composite orbit is a linear 
combination of the actions of the periodic orbits iS„ = ^^g-p^ rripSp. The coefficients rup are 
positive integers and depend on n though we do not indicate this in the notation. More precisely, 
the index n labels all possible finite linear combinations of actions of primitive periodic orbits 
with positive coefficients. In addition there is a zero- length contribution, n = 0, for which 
(So = 0. The zeta function is a sum over all composite orbits in the form 

C{E) = f2Cnexp(^^Sr}j , (4) 



where 



ptv. inSdet(M^-l) 

and Co = 1. Here fip is the Maslov index of the primitive periodic orbit p, and Mp is its stability 
matrix. Mp has eigenvalues apB^^ where cxp is the sign of the eigenvalues. In (jl]) with ([5]) it 
is assumed that all periodic orbits are unstable. We note that in one- dimensional systems the 
corresponding approximation for ({E) consists of a finite number of terms, see [T^ . 

Due to the exponential proliferation of periodic orbits in chaotic systems, the sum in (jlj) 
is not absolutely convergent for real E. Motivated by the Riemann-Siegel formula for the 
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Riemann zeta function, Berry and Keating derived a resummation of the sum based on a 
functional equation of the spectral determinant in h, A{E,h) = A{E,—h) [Ij. (A different 
approach is based on the fact that A{E) is real if E is real ^20].) The derivation uses Cauchy 
integral techniques and the semiclassical sum (jl]) in regions of the complex /i-plane where it 
is absolutely convergent. This resummation makes the sum effectively finite by expressing the 
contribution of longer composite orbits in terms of shorter ones. The approximation has the 
form (in the semiclassical regime) 

AiE) ^ 5(E) 2Ree--^(^) f^C„exp Us^) 9'\Th/2 - %) , (6) 

n=0 ^ ^ 

where 7^ = d5„/ dE are the periods of the composite orbits. 9^^{T) is a function which smoothly 
truncates the sum. It goes from zero to one in a region with width of order h}~^/'^ around T = 
(for the optimal choice of one free parameter). Th is the Heisenberg time Th = 2TThd{E) oc 
h^-^, where d{E) = N'{E) is the mean density of states. In [1] the truncation function is 
given by a series of terms whose leading term contains a complementary error function, and 
the correction terms are given by higher order derivatives of the complementary error function 
and are different from zero only in the vicinity of T = 0. 

Resummed spectral determinants form the basis for the semiclassical approximations in this 
article. The arguments of Berry and Keating apply also to the more general expansions of spec- 
tral determinants that are considered in this article. This is because the spectral determinants 
satisfy the same functional equation, and the semiclassical expansions have the same general 
structure as in ([2]) with (jlj). The difference is that not all of the trajectories that form the 
composite orbits are periodic. 

In the next section we will consider Hamiltonian with singular perturbations. In this case 
there are additional contributions from classical trajectories which are not periodic. They enter 
the semiclassical approximation by means of the Green's function G{q, q', E) = {q\{E~H)^^\q') 
whose semiclassical approximation is given by [fj 

G(g,q',E)^^v^^/^ (7) 

7 

where 

(8) 

Here 7 labels all classical trajectories at energy E from q' to q. is the action and the 
number of conjugate points along 7. The prime at the determinant in ([8]) indicates that the 
determinant involves local coordinates perpendicular to the trajectory. 

The approximation ([7]) is not correct in the limit q' — > q, because it does not describe 
correctly the contribution of those trajectories whose length goes to zero in this limit. The 
correct contribution of these short-length trajectories is given in this limit, for Hamiltonians of 
the form H = /2m + V{q), by the free Green's function [211 122| 

where H^j^_-^ denotes a Hankel function of the first kind [23]. The classical momentum p is 
evaluated at {q + q')/2. Formula shows the divergence of the Green's function as q' ^ q 
for D > 1. The divergence is logarithmic for D = 2 and proportional to l/\q — q'\ for D = 3. 



2n 



(27ri?l){^+i)/2 



^det' ^ 
qq' \dqdq' 
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From the imaginary part of expression ([9]) one can obtain the Thomas-Fermi approximation 
for the local density of states 

oo ^ 

d{q,E) := J2 mq)ME-E^) = -- lmG{q,q, E + ie)^^, , (10) 

771=1 

in the form 

Integration over q and E yields ([3]) for the Hamiltonians considered. 

For convenience in the next sections we introduce the following notation. Let a be any 
quantity which has a semiclassical expansion in terms of classical trajectories or composite 
trajectories in the form a^^ = Yin o-n^^'^"^^- Then we denote the smoothly truncated semiclassical 
sum by 

„sc,tr ^ ^a„e^^n/?^r(TH/2 - t„) , (12) 

n 

where sc stands for semiclassical and tr for truncation. In this notation A(i?) ^ 2 Re A'^'^'*''(ii^), 
and the semiclassical energies satisfy the condition 

2ReA''^'\E'^) = 0. (13) 



3 Singular perturbations 

We discuss in this section singular rank one and rank two perturbations as a means to derive 
semiclassical approximations for wavefunctions and Green's functions in the next section. These 
kind of perturbations can be defined by the applying the method of self-adjoint extension. It is 
instructive to start with the example of a one-dimensional delta-function potential, for which 
the correct results can be obtained by simpler calculations. 
Consider a quantum Hamiltonian of the form 

H = Ho + XA, (14) 

where Hq is the unperturbed Hamiltonian and the perturbation has the form A = \q){q\- The 
eigenstates of H are determined by the condition 

= {Ho-E)\ij) + \\q) {q\ij). (15) 

Multiplying this equation from the left by (g| Go, where Go = {E — Ho)^^ is the resolvent of 
the unperturbed Hamiltonian Ho, yields 

= (l-A(g|Go|g))(g|V^). (16) 

The solutions are given by unperturbed energies E^ if the corresponding eigenstate vanishes 
at q, and by the zeros of the bracket on the right-hand side. 

The resolvent of the perturbed Hamiltonian can be obtained from the general relation 
G = Go + GoXAG from which follows that 

^ = l-Gl\q){q\^' ■ ^'^^ 
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The eigenvalues of the Hamiltonian H are given by the poles of the resolvent G. Hence the 
new spectral determinant has the form 



A{E) = Ao{E) det(l - AGo|g)(g|) • (18) 
Applying detM = e'^'''°s*^ yields 

det(l - \Go\q){q\) = e" ^^-i ^^(^^I'^^^'^D" = g-^^^"-! ^^'^'G"'"))" = 1 - \{q\Go\q) , (19) 
and the spectral determinant follows as 

A{E) = Ao{E)[l-\Go{q,q,E)], (20) 

where Go{q,q', E) = {q\Go\q') is the Green's function. Note that the poles of the Green's 
function cancel the zeros of Ao{E) in (!20l) . if they are non-degenerate and the corresponding 
wavefunction does not vanish at q. The zeros of A{E) are given by the zeros of the expression 
in the bracket, plus the zeros of Aq{E) with a degeneracy that is reduced by one (except if all 
functions in the eigenspace vanish at q). This agrees with the discussion after equation (ITB]) . 
An alternative presentation of the Green's function can be obtained by using the relation 

1+ "l!><l^°,. (21) 



l-GoA|g){,| l-\{q\GM 

which can be verified by multiplying both sides with the denominators. Inserting this into f[T7l) 
results in 

«'^) = «°+«»i«>r^A(kk><'"'°' '''' 

If the fraction is expanded into a geometric series one obtains an expression that can be inter- 
preted as a sum over diffractive orbits [24]. 

It is well known that a delta-function potential is not well-defined in two or three dimensions. 
This is reflected by the fact that (qlGol?) is infinite in these cases. The problem of properly 
defining a pointlike perturbation for D = 2 and D = 3 has a long history, and the extensive 
literature about it is reviewed in [25t 126] . There exist several methods for dealing with this 
problem. The standard method nowadays is to apply the theory of self-adjoint extension to 
a Hamiltonian from whose domain one point has been removed, see e.g. f27\ [251 UE\ [22] • 
The resolvent of the Hamiltonian that is obtained by the self-adjoint extension is related by 
Krein's formula to the resolvent of the unperturbed Hamiltonian. The resulting expression for 
the resolvent of the perturbed Hamiltonian has a very similar form as (!22|) . The difference is 
that the Green's function Go{q, q, E) = {q\Go\q) is replaced by a regularized expression. 

GiE) =Go + Go\q)- r^A-—^Aq\Go , (23) 

and the corresponding expression for the determinant is 

Ai(E, A, q) = Ao{E) [1 - XGl^'^iq, q, E)] . (24) 

The regularization of the Green's function is not unique. Different regularizations can differ by 
an arbitrary real constant [29j. However, all different regularizations lead to the same family 
of self-adjoint extensions of the Hamiltonian. This can be seen from the fact that a change 
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of the regularized Green's function in (123!) by an additive constant c can be compensated by 
redefining the perturbation parameter according to A — > A/(l + cA). 

We will regularize the Green's function in the following by subtracting its divergent part in 
the limit — > q. For Hamiltonians of the form Hq = /2m + V{q) the asymptotic behavior 
of the Green's function in the limit q' — * g is given by the free Green's function ([9]), and we 
define 



q'-*q 



m 



G'^^iq, q, E) = hm G^q, q\ E) - — log(fco|q - q'\) [D 



G'^^iq, q, E) = lim 

q'^q 



G^q,q',E) + 



27rfp\q — q'\ 



:D = 3). (25) 



The constant fco in the first line of fl25p is arbitrary and is included for dimensional reasons. 
It follows from the freedom in the choice of the regularization that the constant A cannot, in 
general, be interpreted as the physical strength of the scatterer. It is interesting to note that a 
possible physical interpretation can be obtained by introducing a new parameter a 

y = -^log(M {D = 2), \ = ^ {D = 3). (26) 
A irn^ A 2Tih'^a 

Then formula (1231) is valid for the scattering from a small but finite scatterer of radius a in the 
limit that a is much smaller than the wavelength |30j . 

The semiclassical approximation for the determinant flMl) is obtained by applying the semi- 
classical resummation that was discussed in the last section. It results in 

Ai(E, \q)^2 M^oiE) (1 - \Gl'^{q, q, E^''^ . (27) 

This has the following interpretation. Ao{E) and (1 — XG^Q^{q, q, E)) both have a semiclassical 
expansion in terms of a sum over classical trajectories. These two sums are multiplied, and one 
obtains a sum over combined trajectories with actions Sn + 5"^. These combined trajectories 
might be viewed as a new form of composite orbits. The corresponding times along the orbits 
are 7^ + T^. The sum over the combined trajectories is then smoothly truncated at half the 
Heisenberg time. 

In the following we will consider also rank two perturbations which correspond to self- 
adjoint extensions of a Hamiltonian from whose domain two points are removed [311 [321 [26] . 
They formally correspond to a perturbation by an operator of the form 

"Ai|qi)(gi| + li\qi){q2\ +Ji\q2){qi\ + Mfl2){q2\' 

where Ai and A2 are real and /i = /i^ + i/ij is complex. In this case the determinant takes the 
form 

A(E) = Ao(^) detF , (28) 
where F is the following 2x2 matrix 



(Gl'^{q,,q,,E) Go{q,,q„E)\ fX, 
\Go{q„q^,E) Gl'^^q^q^E) {fi X, 



(29) 



We will apply the determinant ( [28l) in the following for Ai = A2 = and denote it by 

A^iE, fi, q„ q^) = Ao{E) f{E, fi, q„ q^) , (30) 
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where the function / follows as 

f{E, /i, Qi, Qa) = 1 - f^Go{q2, qi,E)- fiGo{qi, E) 

+ \f^\'[Go{q„ q2, E) Go{q2, ^i, E) - Gl'^q,, q„ E)Gl'^{q2, q^, E)] . (31) 

The determinants A2 and Ai are entire functions of E, because the poles that arise from the 
Green's functions are cancelled by the zeros of Aq. The corresponding semiclassical approx- 
imation for A2 is obtained by inserting the semiclassical expansion for Aq and applying the 
resummation 

A^iE, /i, q^) ^ 2 Re[Ao(E) /(E, /z, q,, q^T''' , (32) 
which is to interpreted according to the discussion after equation fl27|) . 



4 Wavefunctions and Green's functions 

The effect that a perturbation by a pointlike scatterer has on an energy level depends on the 
value of the wavefunction at this point. The semiclassical approach in this section is based on 
this relation. Consider one non-degenerate eigenvalue Em of the unperturbed Hamiltonian Hq. 
The regularized Green's function in the eigenvalue equation f l24l) has a pole a.t E = E^ with 
residue \iprn{q)\'^- Hence we find that the determinant of the perturbed problem, evaluated at 
the unperturbed energy E = E^, satisfies 

A,{Em,X,q) = -A hm AoiE)^-^^^ = -XA',{Em)\^PUq)\' ■ (33) 

We take a derivative with respect to A on both sides of this equation. This leads to an ex- 
pression for the absolute square of the wavefunction as a quotient of derivatives of the spectral 
determinant 



|V'm(q)P 



Inserting the semiclassical approximation fl27j) results in 



^ Ai(E,A,g) 
^Ai(^,A,g) 



(34) 



A=0 

E = Em 



Re[A o{Em) Go{q,q , Em)^''' 
Re[A',iEm)Y 



This is the result for the semiclassical approximation for the modulus square of the wavefunc- 
tion. It gives an approximation in terms of a double sum over composite periodic orbits and 
closed orbits from q to q which is cut off smoothly when the sum of the times along the orbits 
is half the Heisenberg time. The energies E^ in fl35l) are determined semiclassically by f|T3|) . 

Let us discuss the difference to previous semiclassical formulas for the square modulus of the 
wavefunction. Agam and Fishman obtain their approximation by assuming that only closed 
orbits which are almost periodic are important [Hj. They describe these closed orbits in terms 
of nearby periodic orbits and perform a sum over all repetitions of the periodic orbits before 
applying a resummation to the semiclassical sums. Auslaender and Fishman [18] use the fact 
that for tiling billiards on the pseudosphere there exists an exact relation between closed orbits 
and periodic orbits. They do not assume that closed orbits are almost periodic, but perform 
also the sum over repetitions. Equation fl35|) is simpler than the formulas in these articles 
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mainly because of this difference in the resummation. Furthermore, it is not assumed that the 
contributing closed orbits are almost periodic. Other approximations for the square modulus 
of the wavefunction were derived by Poincare surface of section methods [TU ESI UHl [IS]- In 
these cases the resummation is based on the semiclassical unitarity of the transfer operator, or 
the scattering matrix, and is done in terms of the number of intersections with the Poincare 
surface. 

We continue with a further discussion of ( 135|) . One point to note is that ( l35i) is independent 
of the choice of the regularization of the Green's function. This is because the real part of the 
zero-length contribution to the Green's function is multiplied in ( l35l) by Re[Ao(-E'm)]'^'^'*'^ which 
vanishes. 

There is, however, a non-vanishing contribution to (135|1 
zero-length contributions to the Green's function. Using ([10] 
given by 



from the imaginary part of the 
we find that this contribution is 



-vr d(q, Err, 

The derivative of the semiclassical determinant, equation 



(36) 



with 



and 



consists semi- 



classically of two terms (the derivative of B[E) is multiplied by a vanishing term) 



Re[A[,(E^)]^^'*'- = 7rJ(E^) Im[Ao(^^)] 



sc,tr 



-5(E^) Im 



-mN(E.m) 



Tn Cn exp ( -^5, 



Note that d{E) oc h.-^ . 
approximation 



(37) 

If we take only the first term into account we obtain the following 

d{q, Em) 



diEm) 



(38) 



This corresponds to a projection of the classical microcanonical distribution onto the coordinate 
space. When integrated over some small area, this is consistent with the quantum ergodicity 
theorem which, roughly speaking, states that almost all matrix elements of smooth operators 
tend to the microcanonical average in the semiclassical limit [SI [35l [36] . 

Equation ( l35i) can applied for an examination of scars. If a small smoothing of width Ag 
is applied, then the initial and final momenta of the contributing closed orbits differ at most 
by a value of the order of h/Aq, and they can be described in terms of nearby periodic orbits. 
Since the g-dependence in ([35]) is in terms of the semiclassical Green's function, Bogomolny's 
analysis of the structure of wavefunctions near periodic orbits can be applied without additional 
smoothing over an energy window. An effect of the resummation is that the contribution of a 
trajectory 7 is weighted by the semiclassical expression for the spectral determinant, truncated 
at — Ty. The semiclassical criterion for scars of Agam and Fishman [121 [H] can be 

obtained by integrating f[35l) over some tube around the periodic orbits. 

Equation f[35l) does not contain any information about the phase of the wavefunction. It 
is possible to obtain this information by a similar approach. For this purpose we consider the 
determinant ( l30l) for the rank two perturbation. A calculation similar to ( l33l) leads to 



^m(g)V'm(g') = - 



dflr 



+ i 



A2{E,fi,q,q') 



2^A,{E,fi,q,q') 



(39) 



At=0 

E = Em 
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We insert the semiclassical approximation (132|) and obtain 



2 Re [A^(E„ 



V'm(9)^m(g') ~ ^ r.w^ Msc.tr • (^O) 



This is a double sum over composite orbits and trajectories from q' to respectively from q to 
qr', which is truncated smoothly when the sum of the times along orbits and trajectories is half 
the Heisenberg time. Semiclassical approximations for this quantity have been obtained before 
by Poincare surface of section methods, see [15] . (j^Oj) together with (l35l) allows to calculate the 
normalized wavefunction semiclassically everywhere. 

In the case of an degenerate energy level Em equations (1551) and (1101) are formulas for 
{q\Pm\Q) and {q\Pm\(l')i respectively, where is the projection operator onto the correspond- 
ing eigenspace. If the level is (i-fold degenerate one can choose a basis of the eigenspace for 
which d — 1 basis states vanish at q. Hence equations (135|) and (140|) are formulas for the 
remaining state which, in general, does not vanish at q. 

The Green's function can be obtained in a very similar way. From (l30l) follows that 



Ao(i?)Go(g,q',i?)= i 



d . d 
+ 1- 



A2{E,fx,q,q'] 



We insert the semiclassical approximation for the determinant (l32l) and obtain 



(41) 



n ' ^ i^oy-^) '^oy-i^,q,q ) ^ ^oy^) '^oy^^q ,q)\ 
Go{q, q,E)^ „ r a / ^xisc.tr • 



[Ao{E) Go{E, q,q') + Ao{E ) Go{E,q>,q) 
2Re[AoiE)Y 

As in ( HOl) this is a double sum over composite orbits and open trajectories. Resummed expres- 
sions for the Green's function have been obtained by Poincare section methods in [131 HH US], 
and for G^Q^{q, q, E) for tiling billiards on the pseudosphere in [18]. In appendix [B] we discuss 
an example of a one-dimensional system for which formula (H2|) . with a sharp cut-off, yields the 
exact Green's function. 

Given the approximation (H2]) it is straightforward to obtain resummed formulas for other 
quantities that are related to the Green's function. For example, the Wigner transform of the 
Green function is 



Wn 



E) = l d^q' e--^-^'''G,{q + ^q', q - ^q' , E) = h^^l |^ ' (^3) 



where x = {q,p) and Wm{x) is the Wigner function of the m-th eigenstate with energy E^ 
Using ( l42l) and performing the integrals in ( l43i) in stationary phase approximation results in 



2 Re|A„(B)]' 

where the semiclassical approximation to the Wigner function Wq is given by classical trajec- 
tories that satisfy the midpoint rule x = {xf + Xi)/2 where Xf and Xi are the final and initial 
points in phase space [9l [37]. The action terms in (jUl) are of the form Sn + S^ — q'p, where q' is 
determined by the stationary phase condition, but the resummation is done in terms of Tn + T^- 
One may speculate that (Hlj) holds also for a resummation in terms of ^ [Sn + S-y — q'p] , and 
that this can be derived, for example, by starting with a different perturbation. 

The corresponding approximation for the Wigner functions of the eigenstates is given by 



, , [AojEm] Woix,Em) + AoiEm) Wo{x,Em)r''' 
^-^"^ " 2/.-ReK(E„^-- • ^''^ 
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5 Expectation values 



In this section we apply a similar approach to expectation values. Consider the perturbation 
of a Hamiltonian by a self-adjoint operator A in the form H = Hq + XA, and let A{E, A) be 
the associated spectral determinant. As before we denote the determinant of the unperturbed 
Hamiltonian by Ao(-E). In contrast to the previous sections we assume now that A is smooth 
so that the spectral determinant is semiclassically given only in terms of periodic orbits (and 
not diffractive orbits). Then, if -Em (A) is a non-degenerate energy level, and denotes Em{0), 
it follows from the form ([1]) that 



d_ 



A{E,X) 



dE^{X) 



A=0 

E=Err, 



dX 



ME) 



A=0 



E — E„ 



(46) 



where the second equality follows from first order perturbation theory. Hence the expectation 
value of the operator A in the state iV'm) can be written as a quotient of two partial derivatives 
of the spectral determinant 



^A(E,A) 



(47) 



A=0 

E = Err. 



The corresponding semiclassical approximation is given by 

Re[f,A{E,X)r 



{ipmlAllpm) 



Re[^A(i?,A)] 



sc.tr 



(48) 



A=0 

E=Em 



We insert the semiclassical approximation of the determinant ([2]) and perform the derivatives. 
This results in 



{lprn\A\i'm) 



ndAjEj lm[A{Ern)r^'^ - \B{E^) Im [e-'-^(^-) ^^=1 q exp {{S^)] 
7rd{Ej Im[A(E™)]-'f - i5(Em) Im [e--^(^"0 EZiTnC^exp {{Sn)]' 



tr 



Here 



ME) = ^N{E) 



-{2nh)-'> I d^qd''pAiq,p)SiE-Hoiq,p)) 



A=0 



where A{q,p) is the classical function corresponding to the operator A. Furthermore 



Ar, 



dSp 
'dX 



peVn 



A{q{t),p{t))dt 



(49) 
(50) 

(51) 



where the integral is evaluated along periodic orbits. The last equality can be obtained from 
[381I39]. 

Note that both, d^lE) and d{E), are proportional to h^^. If we take only these terms into 
account we obtain the following approximation 



{lpm\A\lprn) 



d{Ej 



(52) 
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This is the average of A{q,p) over the energy shell, in accordance with quantum ergodicity. 

The authors of PD] express expectation values as a quotient of derivatives of the zeta function 
03]) which is modified by including additional terms in the Euler product. For the composite 
orbits these additional terms are factors of the form e"'^^". The resulting approximation differs 
from fH9l) in that zero-length terms were not included, and it does not involve a resummation. 
Instead a curvature expansion was applied. In [33] the scattering approach to quantization 
was applied to obtain expectation values of smooth operators in terms of Poincare section 
wavef unctions. 

6 Discussion 

We have derived in this article semiclassical approximations for wavef unctions. Green's func- 
tions and expectation values in terms of a finite number of classical trajectories by relating 
them to spectral determinants. The main results are equations (135|1 . (HOl) . P2l) . and (148|1 . 
These formulas have a simple and transparent form, and this might facilitate a theoretical 
study of semiclassical properties, for example well numerical evaluation. We have 

concentrated on the position representation, but it is not difficult to obtain other representa- 
tions as well as is shown by the example of phase space representations in equations (jH]) and 
( 1451) . Alternatively, one may apply different perturbations from the start. 

One quantity that appears in almost all approximations as an additional weight of a tra- 
jectory 7 is the semiclassical spectral determinant which is smoothly truncated at half the 
Heisenberg time minus the time along the trajectory, Th/2 — T^. Hence it would be of inter- 
est to understand the dependence of the truncated semiclassical spectral determinant on the 
cut-off. For systems with zero Maslov index one might investigate this dependence by using a 
Guinand-type approach, which involves approximating it by an integral [H]. 

As is demonstrated by the example in appendix [B] the method can also be applied to one- 
dimensional systems. An open question is how the approach of this article is applicable to 
integrable systems with dimension D > 1, or mixed systems. The relations between wavefunc- 
tions, etc., and spectral determinants in equations (IMI) . (15^ . (HTj) . and (H7I) are general and 
do not depend on the nature of the classical dynamics. What is needed is an explicit semiclas- 
sical expansion of the spectral determinant in terms of composite orbits for these systems, or 
a semiclassical approximation of spectral determinants that is based on the EBK-quantization 
for integrable systems. 
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A Regularization of spectral determinants 

This section discusses Weierstrass and zeta regularized spectral determinants and is based on 
[32] • It considers spectral sequences {Ek} with < Eq < Ei < E2 < . . ■ and E^ ^ 00 as 
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— oo. Further assumptions concern the partition function, or the trace of the heat kernel. 
It is assumed that its definition in terms of a sum over energies is convergent for Ret > 0, and 
that it has the following asymptotic expansion 



oo 



e(t):=^e-*^'=~5^Q„f", as t^O, (53) 



k=0 n=0 



with an increasing sequence of real i„ without limit point and io < 0. Let fi = —io in the 
following. Equation (153|) implies that the leading term in Weyl's law for the spectral counting 
function is given by 

There are different ways to define a spectral determinant by a convergent product over the 
eigenvalues E^. of the form ([1]). The Weierstrass canonical product is defined by 

^ f E\ f E E^ ^[^1 \ 

^„,,^).n(l__)e.p(-.^....._^], (55) 

where [fi] is the integer part of fi. An alternative regularization, the zeta regularization, defines 
the spectral determinant by an analytic continuation of the generalized zeta function 

Z(s,a) = ^(Efc + a)-^ = r(s)~M Q{t)e-''H"^dt . (56) 
k=o 

This definition is valid for Res > /i and a > —Eq and is extended elsewhere by analytic 
continuation. The determinant is defined as 

Az{E)=exp[-Z'{0,-E)] , (57) 

where the derivative is with respect to the first variable. Both determinants, Aw{E) and 
Az{E), are entire functions of E. The relation between both determinants can be understood 
by looking at the trace of the /z-th derivative of the resolvent which is given by the convergent 
sum 

oo 

TrG(M)(E) = {-ly^^fAl - Ek^^^"^-' . (58) 

fc=0 

(Note that the resolvent is defined with a different sign in [32]). It can be shown that both 
determinants can be obtained by a multiple integration of 0581) . The difference is that the lower 
integration limit is E = for the Weierstrass regularized determinant whereas it is = — cx) 
for the zeta regularized determinant 

log Aw (E) =U j TrG(M)(^) , logAz(^) =U ) TrG^^H^) • (59) 

In the second case the integral is defined by finite part extractions when necessary. For example, 
the following rules hold 

i-ET dE' = , (s^l), r i-E')-' dE' = - \og{-E) , 



oo 

E 



I log(-E') dE' = E{log{~E) - 1) . (60) 

^— oo 
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It can be further shown that both determinants are related by 



Az{E) = Aw{E)expl-Y,cimEA . (61) 

\ m=0 / 

Expressions for the constants be found in H2] . 

Semiclassical approximations for the spectral determinants are obtained by starting with 
the semiclassical expression for the trace of the [/i]-derivative of the resolvent and doing the 
integrations in flS^ . The zeta regularized determinant is convenient for semiclassical approx- 
imations, because the periodic orbit terms do not give a contribution a.t E = —oo. If the 
periodic orbit terms are integrated for chaotic systems with unstable periodic orbits they yield 
the dynamical zeta function C{E) in (jlj) ^T9[ [3]. The asymptotics of TtG^^^^{E) in the limit 
E — >■ —oo is determined by the zero-length contributions. From the asymptotic expansion of 
the trace of the heat kernel fl53ll follows that 



TrG(M)(E) = -r(H + 1) Z([/i] + 1,-E)r^-J2 c^jm + 1 + ^n){-Er--^^^-' , (62) 

n=0 

as E — s> — oo. For the example of a two-dimensional billiard system with Dirichlet boundary 
conditions (units h = 2m = 1) the first three coefficients are c_i = A/ (47r), c_i/2 = —L/ (Sy^), 
Co = C, and /i = 1. Here A is the area of the billiard, L is its perimeter and C is determined 
by the curvature and corners of the boundary [43j. The integration of (162|) with the rules (l60l) 
then yields the zero- length contributions to the spectral determinant. We state here the result 
for two-dimensional billiard systems 

Az{E) ^ exp {AE{\og{E) - 1)/Att + C\og{E) - i7TN{E)}C{E) , (63) 

where 

N(E) ^ —E - —y/E + C . (64) 

^ ^ 477 477 ^ ^ 

Equation fl63l) is in agreement with the appendix of [7j. Comparison with ([2]) shows that in the 
approximaton fl63|) for the zeta regularized determinant the function B{E) is given by 

Bz{E) = exp{AE{\og{E) -l)/4n + C\og{E)} . (65) 

It is part of the zero- length contribution. For the Weierstrass regularization B{E) contains an 
additional factor according to ( !6T|) . 



B A one-dimensional example 

In this section we apply the resummation method to a one-dimensional example. In order 
to test the formula for complex-valued Green's functions we consider a system without time- 
reversal symmetry. It will be seen that the resummation is exact if we choose a sharp cut-off 
by a Heaviside theta function. 

fO if T<0 

e''{T) = 9{T) = h if T = (66) 
[l if T>0. 
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The Schrodinger equation is given in units where h = 2m = 1 by 



(67) 



where a is constant, and we apply periodic boundary conditions of the form iplq + I) = ip{q). 
The solutions for ( |671) are given by 



n e Z. 



(68) 



The corresponding free Green's function is obtained from its differential equation [E = k"^) 

1 



(69) 



The Green's function of the considered system follows by the method of images 

oo 

Go{q, q', E)= + <l'^ E) . (70) 



This representation of the Green's function is identical to its semiclassical approximation. The 
sum can be interpreted as sum over all trajectories from q' to q with lengths L„ = \q + nl — q'\. 
The sum in flTOj) involves only geometric series and can be carried out with the result 



^ , ^ e-'"^''-^') rsin(a/)e''=(«"«') + sin(A;/-A;|g-g'|)-sin(a/-A;|g-g'|)l 
Go{q,q,E) = ^ — — — — (71) 



2k (cos(aZ) — cos{kl)) 



For q' = q we obtain 



1 

Goiq, q^E) = — + Y, T^e'^"' cos(an/) 

n=l 



sin(A;/) 



2k (cos(a/) — cos(A;/)) 



(72) 



From the zero-length term we obtain the mean spectral staircase N{E) = kl/ir. Instead of the 
Heisenberg time it is more convenient to speak of the Heisenberg length which is Lh = 2kTH = 
AkirN'^E) = 21. Hence the cut-off for the orbit sums is at /. 
The spectral determinant is 



Ao(^) = exp / dETrGo(^) = cos(a/) - cos{kl) . 



This can be written in the form 



(73) 



(74) 



which corresponds to the sum over composite orbits in this one- dimensional system. Applying 
the resummation with the theta-function as cut-off we find that 



Ao(^) ^ 2Re[Ao(^)]'''='*'^ = 2Ree 



-\kl 



i + 1 cos(«Oe'^' 



cos(a/) — cos{kl) . (75) 
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Hence the resummation is exact. 

Let us apply now the resummation to the Green's function 



r ( . p^ ^ [Ao(i^) Go(g, g^ E) + ^,{E) Gojq', g, E)^^'^ 

^° ^ ' 2Re[Ao(i?)]-'- • ^^^^ 

One has to multiply the sums (ITOll and (17^ and discard all terms for which the joint lengths 
of the trajectories is larger than Lh/'^ = I (and multiply terms for which the joint lengths is I 
by 1/2). We consider first the case q' = q and obtain 



Go{q,q,E) = 2Re 



_i_e-i^' [-\ + 1 cos(a/)ei'=' - ie*'^'-*"' - l^e^i+iaii^ 



cos(a/) — cos{kl) 

sin (A;/) 



2k[cos{al) — cos(A;/)] 
which agrees with (1721) . In the case q' q 



(77) 



Go{q,q',E) = 
+ 



J_ -ikl \ _Qik\q~q'\-ia(q-q') _ ^ik{l-\q-q'\)-ia(q-q'Tl)] 
4ifc L J_ 

cos(a/) — cos{kl) 

Lpifc' \ _Q-ik\q-q'\-ia(q-q') _ Q-ik{l-\q-q'\)-ia{q-q'^l)~\ 

ilk L J 



cos(q;/) — cos{kl) 

^ sin(fc/ -k\q-q'\) + sm{k\q - g^|)e^'"' ^-ia(q-q') /^gx 
2A;[cos(a/) - cos(A;/)] ' ^ ^ 

where the upper and lower signs in the exponents correspond to the cases q > q' and q < q', 
respectively. A short calculation shows that (1751) is identical to (17T1) . 
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